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ABSTRACT 



We present high resolution imaging of the nucleus of NGC4258 from 1 /im 
>■ ■ to 18 /xm. Our observations reveal that the previously discovered compact 

{>^ ■ source of emission is unresolved even at the near- infrared resolution of ~0.2" 

5;0 ' FWHM which corresponds to about 7 pc at the distance of the galaxy. This is 

consistent with the source of emission being the region in the neighborhood of 
^ ■ the purported 3.5x10^ Mq black hole. After correcting for about 18 mags of 

visual extinction, the infrared data are consistent with a F,y oc //-i-4±o.i spectrum 
CL|' from 1.1 /im to 18 /im, implying a non-thermal origin. Based on this spectrum, 

O ' the total extinction corrected infrared luminosity (1—20 /im) of the central 

-t— » ' source is 2x10* Lq. We argue that the infrared spectrum and luminosity of the 

ci ■ central source obviates the need for a substantial contribution from a standard, 

thin accretion disk at these wavelengths and calculate the accretion rate through 
^ ' an advection dominated accretion flow to be M ~ 10"'^ Mq/jt. The agreement 

j^ ■ between these observations and the theoretical spectral energy distribution for 

advection dominated flows provides evidence for the existence of an advection 
dominated flow in this low luminosity AGN. 
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1. Introduction 

NGC4258 (M106) is a weakly active Seyfert 1 (Filippenko 1996) galaxy at a distance 
of ~7 Mpc (we adopt Ho=65 km s~^ Mpc~^; at this distance 1"=33.9 pc). The nucleus of 
the galaxy is a prototypical low luminosity AGN. Dynamical evidence for the existence of 
a massive black hole in the nucleus comes from high spatial resolution VLBI observations 
of water masers, which reveal the existence of a thin edge-on disk of inner radius 4 
milliarcseconds and outer radius 8 milliarcseconds (Miyoshi et al. 1995). The high velocity 
features of the masers were found to be symmetrically offset from the systemic velocity by 
700—1000 km s~^ as would be expected in a Keplerian disk. The inferred binding mass of 
3.5±0.1xl0^ Mq within a radius of 0.14 pc (Herrnstein, Greenhill, & Moran 1996) makes 
it the strongest supermassive black hole candidate after the nucleus of the Galaxy (Genzel 
et al. 1997, Ghez et al. 1998). The adopted binding mass and distance for the galaxy are 
consistent with the more recent measurements of Herrnstein et al. (1999). Observations 
by the X-ray satellite ASCA indicate the existence of a strongly obscured source with a 
photon index of r=1.8±0.3 at the nucleus of the galaxy (Makishima et al. 1994). The 
photon index r=H-7 where 7 is the spectral index. The line of sight column density of 
Nj:/=1.5xl0^^ cm~^ derived from a fit to the absorbed X-ray data, corresponds to Ay~80 
mags which would obscure the nucleus at visible wavelengths. The extinction-corrected 
X-ray luminosity from 2—10 keV is 4x10^° ergs s~^ and is relatively insensitive to the 
absorption column density. Based on arguments that the X-ray luminosity is typically 
about 5% to 40% of the bolometric luminosity in AGN (Mushotzky, Done, & Pounds 1993), 
LboI for the central engine in NGC4258 is of order 10^^ ergs s^^. This is ~2xl0^^ Le, where 
Le is the Eddington luminosity of the massive black hole at the nucleus and corresponds to 
4.4 xlO''^ ergs s"^ 

The low bolometric luminosity of the central engine can be explained by an incorrect 
estimation of Lboi since the overall spectral energy distribution of the central engine is not 
well known. Alternatively, since \jboI=ii^(^ where 77 is the radiative efficiency and M is the 
mass accretion rate of the central engine, a low bolometric luminosity could be attributed 
to: 

1. a sub-Eddington mass accretion rate through a standard, geometrically thin disk. For 
the black hole in NGC4258, the Eddington mass accretion rate M£;~0.75 MQ/yr where 77 
has a fiducial value of 10%. 

2. a low radiative efficiency for the accretion phenomenon whereby most of the viscously 
generated energy is carried through the event horizon of the black hole. This is more 
commonly called an advection dominated accretion flow (ADAF; Narayan & Yi 1995). 

Strong line of sight extinction to the nucleus of NGC4258, makes estimating the 
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intrinsic spectral energy distribution difficult. Wilkes et al. (1995; hereafter W95) detected 
a polarized continuum with F^ oc //^i-i±o.2 ^^ scattered light at visible wavelengths. The 
luminosity of the central source is only weakly constrained by these observations to lie in the 
range 10^^ — 10'*^ ergs s~^ and is dependent on the scattering model. Infrared observations 
are less susceptible to extinction effects and are therefore a better probe of the energy 
spectrum and the output luminosity. Chary & Becklin (1997, hereafter CB97) detected a 
compact infrared source in the nucleus of the galaxy. Based on the visible light polarimetric 
observations of W95, CB97 concluded that the intrinsic spectrum of the infrared source is a 
power law (Fj, oc u^^) and that the observed J—H, //—iT colors could be explained by about 
17 mags of visual extinction along the line of sight. While the absolute extinction corrected 
luminosity at these wavelengths is consistent with an ADAF model (Lasota et al. 1996), 
the derived spectrum poses some problems because at these wavelengths, both the ADAF 
model and the standard thin disk model predict a z/3 Comptonized blackbody spectrum 
arising from a cool geometrically thin, optically thick disk (Pringle 1981). In addition, 
infrared emission from thermal reprocessing of the central visible/UV continuum by dust 
in a flared disk or 'torus' could not be unambiguously ruled out. The only observation at 
mid-infrared wavelengths which can potentially trace the dust distribution or constrain the 
overall spectrum of the central source is the 8—13 /im 5.7" diameter data from Rieke & 
Lebofsky (1978; hereafter RL78). This measurement could be contaminated by extended 
emission from regions of nuclear star formation. 

Lasota et al. (1996) applied an ADAF model to fit the X-ray spectrum of NGC4258. 
The model consists of an inner, optically thin fiow which is geometrically thick. The 
ADAF exists from an inner radius of 3 R5 up to a radius of 10—100 R5 where R^ is the 
Schwarzschild radius of the central black hole and corresponds to 3.4x10^^ pc. The ions 
which are heated by the viscously generated energy have a temperature of ~10^^ K but 
their radiative timescale is longer than the accretion timescale resulting in most of the 
energy being transported into the event horizon of the black hole by radial convection or 
'advection'. The electrons (Te~10^ K) produce synchrotron photons which could be the 
seeds for self-Comptonization to produce a hard X-ray spectrum. The model also includes 
a cooler geometrically thin, optically thick disk which extends up to 10^ Rg to reproduce 
the blackbody bump seen at visible/UV wavelengths in most AGN. This kind of accretion 
flow allows a 'normal' accretion rate of m~0.2a but requires a low radiative efficiency since 
most of the viscously generated energy is transported into the event horizon of the black 
hole by the advection process. Here, a < 1 is the dimensionless viscosity parameter defined 
by Shakura & Sunyaev (1973), while rh = M{MQ/yr)/ Me is the accretion rate in units 
of the Eddington rate. Mahadevan (1997; hereafter M97) has investigated the spectral 
energy distribution from an ADAF in detail and shown that the Comptonization of soft 
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cyclosynchrotron photons results in a power law spectrum between 10^^ and 10^'^ Hz. The 
spectral index of the power law is essentially constrained by the mass accretion rate through 
the ADAF. 

In this paper, we present high resolution imaging of the nucleus of NGC4258 at 
mid-infrared and near-infrared wavelengths. The observations place a better handle on the 
spectrum and bolometric luminosity of the central engine. This allows us to constrain the 
mass accretion rate M of the black hole based on the different accretion models which have 
been proposed. 



2. Observations and Results 

The 1.1 /xm to 2.2 yum data were obtained using the Near Infrared Camera and Multi- 
Object Spectrometer (NICMOS) Camera 2 on the Hubble Space Telescope (Thompson et 
al. 1998). The camera has a 256x256 HgCdTe array with pixel scales of 0.0762" and 0.0755" 
per pixel in x and y, respectively, providing a ~ 19.5" x 19.3" field of view. Frames were 
taken in three filters: FllOW, F160W and F222M, which correspond to 1.1 /xm, 1.6 /im, and 
2.2 /xm respectively. Observations of the galaxy were made using a four-point spiral dither 
in each filter setting; the offset per dither position was 25.5 pixels (1.9"). Non-destructive 
reads (MULTIACCUM) were obtained at each dither position, with integration times of 
40 (FllOW, F160W) and 56 (F222M) seconds per position. The total integration times 
were thus 160 (FllOW, F160W) and 224 (F222M) seconds. In addition to the galaxy 
observations, blank-sky observations were done with the F222M filter to determine the 
background level at 2.2 /im, and dark frames were taken using the same MULTIACCUM 
sequences as the galaxy observations. 

The data were reduced with IRAF. The dark was first created, then the NICMOS data 
were dark subtracted, fiatfielded and corrected for cosmic rays using the IRAF pipeline 
reduction routine CALNICA (Bushouse 1997; See Scoville et al. 1999 for details about 
the reduction procedure). The calibrated images contained pixels with reduced quantum 
efficiency due to contaminants on the array surface, thus a mask was created to minimized 
their effect. The coronographic hole and column 128 were also masked on all images, the 
latter feature was masked due to its sensitivity to minute discrepancies in dark subtraction 
(i.e., the dark current rises sharply toward the center column of the array). The dithered 
images were then shifted and averaged using the DRIZZLE routine in IRAF (Fruchter & 
Hook 1997). Finally, a background level of 0.45 ADU/s was measured from the blank-sky 
data, and this amount was subtracted from the 2.2 /xm image of NGC 4258. The background 
level at 1.1 /im and 1.6 /im was ignored since it is typically less than 0.01 ADU/s. The 
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plate scales of the final 'drizzled' images, shown in Figure 1, are 0.0381" and 0.0378" per 
pixel in x and y. 

Flux calibration of the images were done using scaling factors of 2.195x10^^, 
2.207x10"^ and 5.583x10"^ Jy (ADU/sec)-i at 1.1, 1.6, and 2.2 /xm (Rieke et al. 1999). 
The corresponding magnitudes were calculated using the zero-points 1909, 1087 and 665 Jy 
at mag for 1.1, 1.6 and 2.2 yum respectively. The 1.1 /im, 1.6 ^m. and 2.2 /xm images were 
aligned to within 0.2 pixel using centroids and elliptical isophotes generated. Synthetic 
NICMOS point spread functions (PSFs) were generated using the TINYTIM routine (Krist 
& Hook 1997) since no observed PSFs were available. 

As shown by CB97, the nucleus of the galaxy is clearly redder than the surrounding 
starlight. The colors of the nuclear regions (~0.12" radius from the center) are 
[1.6] — [2.2]=0.9±0.1 mag and [1.1] — [1.6]=1.1±0.1 mag. In comparison, the colors of 
the galaxy in an annulus of radius 0.5" to 1.0" are [1.6] — [2.2]=0.52±0.05 mag and 
[1.1] — [1.6]=0.91±0.05 mag. The reddening can be attributed to either a red central source 
or dust extinction of starlight. CB97 argued for the former and we provide additional 
evidence in this paper that this is true. To estimate the contribution of the putative 
source seen by CB97 at the nucleus of the galaxy, it was necessary to subtract a stellar 
background from the galaxy frames. The 1.1 /im surface brightness profile was fitted with 
an exponential and subtracted off such that the residual consisted of a central core with 
a FWHM corresponding to the FWHM of the 1.1 /im PSF (~0.1"). There is no evidence 
of an Airy ring in the residual and hence we conclude that we do not detect a compact 
component of emission at 1.1 /xm. The residual has a flux density of 0.5 mJy and we 
consider it as an upper limit to a point source. 

Since the 1.1 /xm data do not show a compact source, it is reasonable to assume that 
it traces the distribution of starlight in the central regions of the galaxy. Furthermore, the 
contribution from the polarized component observed by W95, based on it's flux density and 
spectrum, is found to be less than 0.1 mJy at 1.1 /im. A scaled 1.1 /im surface brightness 
image was subtracted from the 1.6 /im and 2.2 /im surface brightness images. The scale 
factor was derived by performing an intensity-weighted least squares fit within the central 
arcsecond of the form: S[2.2]=ai ^[i.i] + bi P[2.2] where S is the surface brightness, ai 
and bi are constants and P is the TINYTIM generated PSF. For the 2.2 /im image, we 
obtain ai=1.25±0.05 and for the 1.6 /xm image, we get ai=1.29±0.05. This corresponds to 
[1.1] — [2.2] and [1.1] — [1.6] colors of 1.39±0.05 mag and 0.89±0.05 mag respectively. 

The Airy ring is clearly visible in the starlight subtracted 1.6 /im and 2.2 /im images 
which were obtained as described above. The FWHM of the residual is also very similar to 
the FWHM of the synthetic PSFs at the respective wavelengths. We thereby conclude that 
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the compact component of emission at 1.6 ;um and 2.2 /im is unresolved at our resolution 
of ~0.2". This implies that the size of the emitting region has a radius ;^3.5 pc at the 
distance of the galaxy and is consistent with the emission originating from the vicinity of 
the black hole at the nucleus. The profile of the central source is shown in Figure 2. The 
brightness of the central source was measured by standard circular aperture photometry. 
The observed flux density of the compact source is 4.0 ±0.7 mJy at 2.2 /xm and 0.9 ±0.3 
mJy at 1.6 ;um. This agrees well with the values of 4.5 mJy and 1.1 mJy obtained by CB97 
at K and H respectively. 

Rieke (1999) has suggested that using the standard NICMOS calibration numbers 
for red objects as above, could result in an error of as much as 0.7 mags in the intrinsic 
NICMOS magnitude especially at [1.1]. The effect is smaller at [1.6] and negligible at 
[2.2]. This is because red objects would shift the effective wavelength and the zero point 
of the FllOW and F160W filters, both of which have passbands much wider than the 
ground-based J and H filters. We estimate the amount of this error for NGC4258 by 
applying a conversion from instrumental NICMOS magnitudes to ground-based magnitudes 
(Appendix A). This also serves a secondary purpose which is to compare the colors of the 
starlight in NGC4258, derived from the NICMOS data, with ground-based colors. 

We have converted the NICMOS [1.1] and [1.6] images to ground-based J and H 
magnitudes and rederived the compact source fiux density using the procedure described 
earlier. We find that the starlight in the galaxy has a J— /^ color of 1.0 mag and a J—Hco\oi 
of 0.6 mag. This is in excellent agreement with the ground-based colors of J—K=Q.9 mag 
and J—H=0.6 mag obtained by CB97 at a distance of 6" from the center. These colors are 
consistent with the colors of late type stars expected to dominate the surface brightness 
in the central regions of the galaxy. The resultant fiux density of the central source is 1.1 
mJy at 1.6 fim. and 3.6 mJy at 2.2 /im. This is within the photometric errors assigned to 
the compact component above and hence we conclude that the wide bandpass of the filters 
doesn't significantly affect our determination of the compact source flux density. 

Near infrared data were also obtained in the Lg (Aq = 3.45 /im, AA = 0.57 /im) and 
Kg (Ao = 2.15 /im, AA = 0.33 /im) bands using the Cassegrain IR Camera on the 5.1 m 
Hale telescope. A Kg PSF source was scaled to the flux density of the compact 2.2 /xm 
source in NGC4258 and subtracted from the Kg image of the galaxy after aligning. The 
residual image corresponds to the surface brightness proflle of the starlight in the galaxy. 
This was then normalized to the Lg surface brightness at a radius of 2.25" and subtracted. 
The normalization factor of 'Eks/^Ls=^-Q8 corresponds to a Kg—Lg color of 0.16 mag, 
which is reasonable considering the J— i/and iJ— if colors of the starlight. Photometry was 
performed on the Lg residual image which is dominated by the compact component. We 



measure a flux density of 20±3 mJy for this source. 

The mid infrared data were obtained using the MIRLIN camera on the 10 meter Keck 
II telescope in March 1998 (Ressler et al. 1994). The array is a 128x128 array with a pixel 
scale of ~0.14"/pixel at the f/40 bent Cassegrain focus of Keck II. The instrument has a 
square wave chopper which chops at 5 Hz. The chopper throw was set at 6" which allowed 
the source to be on the array at both chop positions. Observations were made with a filter 
which is centered at 12.5 /im with a passband of 1.2 /im and a filter which is centered at 
17.9 lira with a passband of 2 /im. We chopped 300 times per position and nodded by 6" 
after a minute. The total integration time was 300 s at 12.5 /xm and 450 s at 17.9 /im. The 
resultant la sensitivity at 12.5 //m and 17.9 /im was 20 mJy and 50 mJy in a 2.8" diameter 
beam respectively. 

The mid infrared data were reduced in a standard manner (Bock et al. 1998). Airmass 
corrections were performed which were less than 5% since the standards and source were 
at similar elevation angles. The final image at these wavelengths is dominated by a single 
unresolved source (FWIIM^O.5") with no evidence of extended emission on scales of 
diameter up to 6". The observed flux density of the unresolved central source is 165±20 
mJy at 12.5 /im and 300±30 mJy at 17.9 /im. 



3. Emission from a Compact Extincted Continuum Source 

CB97 showed that the red near- infrared colors of the nucleus of the galaxy could not 
be explained by uniform foreground extinction of starlight or by stars mixed with dust. In 
addition, starlight cannot explain the observed spectrum of the enhanced nuclear emission 
at A > 3 /im (Table 1). 

If the enhanced infrared emission in the nucleus arises in a compact nuclear continuum 
source associated with the 3.5x10^ Mq black hole as argued by CB97, it implies that the 
extinction observed in the X-ray data occurs in a foreground screen such as a dust lane, 
attenuating both the near and mid infrared intensities^. The visible/UV polarimetry data 
on the central source shows a Fj, oc //-ii^o-^ spectrum in scattered light (W95) similar 
to our observed mid-infrared spectrum. Therefore, it is reasonable to assume that the 
intrinsic near-infrared spectrum of the compact source follows a similar trend and that 
the observed [1.6] — [2.2] color of the compact source is due to extinction. We calculate 



^For the rest of this paper, we assume A[ii]=0.282 Ay, A[ig]=0.624 A[ii], A[2.2]=0.382 A[ii], 
A[3.5]=0.182 A[i.i], A[io]-0.192 A[i.i], A[i2.5]=0.098 A[i.i], A[i7.9]=0.083 A[i.i] (Mathis 1990). 
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that the amount of extinction required to fit the near and mid infrared data by a single 
power law corresponds to Ay =18 mag. This is smaller than the extinction derived from 
the X-ray column density. However, as pointed out by Gammie, Narayan, & Blandford 
(1999), the column density Nh derived by Makishima et al. (1994) is unreliable since it 
is sensitive to the relative contributions of the components used to fit the X-ray data. In 
addition, the conversions from the X-ray opacity, tx to Nh and from Nh to Ay is not well 
established. The resultant extinction corrected infrared spectrum follows a z/^i-4±o.i power 
law, suggesting a non-thermal origin. The extinction corrected fiux densities are shown in 
Figure 3 and Table 1. 

The extrapolated fiux density of the central source at 1.1 fim. based on our fitted 
spectrum is calculated to be 10 mJy. If we extinct this 1.1 /im source by Ay=18 mag, we 
obtain a fiux density of 0.1 mJy, consistent with our derived upper limit. The 0.6—20 /im 
luminosity of the central source derived from our fitted F^ oc z/~^'^ spectrum is 8.5x10'^^ 
ergs s~^. This is a factor of 20 greater than the X-ray luminosity. 

RL78 measured a fiux density of 100±12 mJy in a 5.7" diameter beam with a 8—13 /xm 
filter. Since our mid-infrared images do not show any evidence of extended emission 
from the nucleus of the galaxy at scales of 6", this measurement must be dominated by 
the central source. If our extinction corrected measurements are interpolated to 10 /im 
assuming a z/"^'^ spectrum, we obtain an unextincted fiux density in the 8—13 /xm filter 
of 190 mJy. After putting in the extinction in the 8—13 /im window including the 9.7 /xm 
silicate feature, corresponding to Ay=18 mag, we derive an extincted 10 /xm fiux density of 
98 mJy for the central source, in excellent agreement with the large beam measurements. 

The VLA observations of Turner & Ho (1994) had detected the existence of a 
non-thermal source in the nucleus of the galaxy at 2 cm and 6 cm. Milli-arcsecond 
resolution VLBI continuum observations indicate that most of the non-thermal radio 
emission arises in a sub-parsec-scale jet which is aligned with the rotation axis of the masing 
disk (Herrnstein et al. 1997). The black hole region itself is undetected and Herrnstein 
et al. (1998) obtain a 3a upper limit of 220 /iJy at 22 GHz. When compared to the 
extinction corrected infrared spectrum, this seems to suggest that the spectrum has a 
turnover probably due to synchrotron self-absorption between 20 /im and 0.2 mm. The 
total luminosity of the central source longward of 1 /xm is relatively insensitive to the 
turnover frequency and is constrained to lie between 8x10^^ ergs s~^ and 3x10*^^ ergs s~^ 
for turnover wavelengths of 20 /xm and 0.2 mm respectively. Spinoglio & Malkan (1989) 
suggested that the 12 /xm luminosity is typically about 1/5 of the bolometric luminosity of 
active galaxies, independent of whether the emission is thermal or non-thermal in origin. 
We find by this argument, the bolometric luminosity of NGC4258 to be ~10^^ ergs s~^. 
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which is in between our independently derived luminosity range. 

As indicated by the extinction correction to our measurements, the nucleus of NGC4258 
clearly has dust associated with it. The extinction corrected infrared emission from 1 fiui to 
20 fira has a z/^^*^ power law spectrum which appears to extend up to visible wavelengths 
as seen in the polarimetric observations of W95. This makes thermal dust emission at these 
wavelengths unlikely and we infer the infrared emission is non-thermal in origin. Many 
previous authors have concluded that the source of both the near and mid infrared emission 
in the nucleus of Seyferts is thermal emission from dust reprocessing of a blue, nuclear 
continuum (e.g. Giuricin, Mardirossian, & Mezzetti 1995). This was based on a weak 
correlation between the X-ray and mid-infrared emission in a large sample of Seyferts. Our 
derived spectrum suggests that it is more likely that the dust in the nucleus of NGC4258 is 
cool and contributes substantially only at far-infrared/submillimeter wavelengths. 



4. The Accretion Flow in NGC4258 

4.1. Implications for an Advection Dominated Flow 

M97 has derived the 10® to 10^° Hz spectrum of systems dominated by an ADAF 
(Figure 4). The spectrum consists of three components; a synchrotron component at 
u < 10^^ Hz, a power law component for 10^^ < u < 10^^ Hz from the inverse Compton 
scattering of the cyclosynchrotron photons and a bremsstrahlung component at higher 
frequencies. 

From M97, the spectral index of the power law (L,^ oc z/"''') is given by 



7 = 



In A 



„ , /0.3\ /0.5\ / 3 \^/2 
(23.87 m) — — (2) 



A = 1 + 4^e + 16^e (3) 

where Tes is the optical depth to electron scattering, A is the mean photon energy 
amplification factor in a single scattering, ^g = {kTe) / {rrieC^) , C\ = 0.5 is a term relating the 
radial velocity of the accretion flow to the Keplerian velocity and Tmin is the inner radius of 
the ADAF which we assume to be 3 R5. Te is the electron temperature which is constant 
over the ADAF region. Based on the figures from M97, we adopt a Te=4.4xl0^ K, q;=0.3, 
and derive from the above equations, an accretion rate of m ~ 1.2 x 10~^ in Eddington 
units, which corresponds to M ~ 9 x 10""^ Mq/jt. 
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lt is also possible to estimate m from the bolometric luminosity as shown by M97. 
For a=0.3, (3=0.5, we find that our LboI ~ 10^^ ergs s^^ corresponds to m ~ 6 x 10^^ in 
Eddington units which is again reasonable considering the factor of 3 uncertainty in the 
bolometric luminosity. (3 is the ratio of gas pressure to total pressure. 

In comparison, Lasota et al. (1996) derived a rate of m ~ 0.16a through an ADAF. 
However, at that time the infrared data which better constrain the bolometric luminosity 
of the central engine was not available and the ADAF models were fit to the slope of the 
X-ray data. More recently, Gammie et al. (1999) derive an accretion rate of 0.01 Mq/jt 
through an ADAF using the infrared data of CB97 as a constraint. Again, this analysis 
assumed that most of the visible/IR emission has a u^^^ Comptonized blackbody spectrum 
and arises in the cool, thin disk rather than in the ADAF itself. Our recent data suggests 
that the visible/IR emission has a common origin and has an intrinsic spectrum that can 
be well represented by a ~z/^^ power law. 

Herrnstein et al. (1998) have shown based on their 22 GHz upper limit, that the 
ADAF cannot extend significantly beyond ~10^ R5. We find that the 22 GHz limit in 
conjunction with the spectrum of M97 in the self-absorbed radio regime is inconsistent 
with our measurements (Figure 4). We predict that the radio synchrotron spectrum must 
have a spectral index of 2.5> 7 > 1.2 where F^ oc u^ rather than the 7 = 1 derived 
for a thermal distribution of electrons by M97. This could arise if the distribution of 
electrons has a power law tail at higher energies. Maintaining a power law distribution 
of relativistic particles without resorting to an independent injection mechanism such as 
magnetic reconnection has been an unsolved problem for some time now (McCray 1969). All 
previous work has assumed that the particle distribution in an ADAF is thermal. However, 
Mahadevan & Quataert (1997) investigated the emission spectrum for different accretion 
rates through an ADAF. They found that for low accretion rates, the thermalization of 
the electron distribution through Coulomb collisions with the protons and through self 
absorption of synchrotron photons, is inefficient. Adiabatic compression of the electron gas 
acts as the primary heating mechanism which results in the high energy tail of the electron 
distribution function having a non-thermal Gaussian profile rather than a Maxwellian 
profile. Alternatively, Mahadevan et al. (1997) has suggested that if the weak Galactic 
center source Sgr A* is associated with the EGRET source 2EG J1746-2852, the gamma ray 
spectrum could be explained by an E~^'^^ distribution of protons in an ADAF. A similar 
power law distribution of electrons would indicate that the electron distribution function 
is generated by the same viscous heating mechanism responsible for generating the proton 
distribution function. 
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4.2. Implications for a Standard Thin Disk 

We find that the excellent agreement between the fit F^, oc z/~^'^ spectrum and data 
obviates the need for a significant z/3 Comptonized blackbody continuum from a standard 
optically thick, geometrically thin disk. However, 22 GHz observations show the existence 
of a thin disk with inner radius 0.14 pc which sustains the water masers (Miyoshi et al. 
1995). From standard disk theory (Frank, King & Raine 1992), 

T,ff ~ 3.5 X 10^ m-i/^ m^/^ r'^/^ K (4) 

where m is the central mass in units of Mq and r is the radius of the disk in units of 
R5. For a geometrically thin, optically thick disk to be capable of emitting blackbody 
near-infrared photons, it should be at least 1500 K. Based on the derived accretion rate, 
this implies a conservative lower limit of 200 R5 for the inner radius of the thin disk. 



5. Conclusions 

NGC4258 has associated with it a compact nuclear source at near and mid infrared 
wavelengths which is unresolved even at a near-infrared resolution of ~0.2" FWHM. This 
implies the emitting region is ^3.5 pc in radius at the distance of the galaxy. The observed 
infrared flux densities can be flt by a single power law of the form Fj, oc //-i-4±oi from 
1 /im to 18 /im after correcting for about 18 mags of visual extinction. The total extinction 
corrected infrared luminosity (1—20 fj,m) of the central source is 8x10^^ ergs s^^. At shorter 
wavelengths, the power law spectrum seems to extend upto visible wavelengths as shown 
by the polarimetric observations of W95. At longer wavelengths, the infrared spectrum 
along with the 22 GHz upper limit points towards a synchrotron self-absorption turnover 
wavelength between 20 /xm and 0.2 mm. From this spectrum, we derive the bolometric 
luminosity of the nuclear source to be ~2xl0^^ ergs s~^. The low bolometric luminosity 
and power law spectrum extending from visible to mid-infrared wavelengths indicates 
that self-Comptonized synchrotron emission arising in an ADAF is responsible for the 
enhanced infrared emission. The spectrum obviates the need for a "big blue bump" at 
these wavelengths from a cool, optically thick, geometrically thin disk extending beyond 
the advection dominated flow. Such a disk is predicted to have a Fj^ oc u^'^ spectrum which 
appears unlikely from our data. Based on the spectrum of the source and the derived 
bolometric luminosity, we infer the accretion rate through an ADAF in NGC4258 to be 

I©/ 



40-3 M0/yr. 
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A. Conversion of NICMOS magnitudes to Ground-based magnitudes 

In this section, we describe the color correction applied to convert the NICMOS [1.1], 
[1.6] and [2.2] magnitudes to ground-based J, H, K magnitudes. The color correction is 
applied for two reasons: 

1. To compare the results from this NICMOS data set with previous ground-based 
observations, 

2. To correct for calibration errors which arise due to the wider bandpass of the NICMOS 
FllOW and F160W filters. It has been suggested that since the NICMOS filters were 
calibrated off solar analog stars, using the standard NICMOS calibration numbers could 
result in an error of as much as 0.7 mags in the intrinsic NICMOS magnitude for red 

[J — H > 1 mag) objects (Rieke 1999). This error is not significant in the ground-based J, 
H filters because of their much narrower bandpass. 

NICMOS has observed a sample of standards of different colors for which the ground- 
based magnitudes are known (Table 2). Thus, by converting instrumental magnitudes to 
ground-based magnitudes on the basis of these observations, we can eliminate calibration 
errors. Figure 5 shows the color correction with the best corresponding linear fit listed 
below. The color correction at [2.2] is small enough that we take it to be 0. 

J -[1.1] = -0.58 X ([1.1] - [1.6]) + 0.11 (Al) 

H-[l.6] = -0.24 X ([1.6] - [2.2]) + 0.01 (A2) 
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20" 



Fig. 1. — Clockwise from top left: Morphology of the nuclear regions in NGC4258 at 
1.1 /im, 1.6 /im, 12.5 /im and 2.2 /im. The near-infrared images are 20" on a side. 
The 1.1 /im and 1.6 /im images clearly show the dust lane in the nuclear regions of the 
galaxy. The 12.5 /im image which is 6" on a side does not show any extended emission. 
a = 12h 18m 57.5s, 5 = 47° 18' 14" (J2000). 
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Fig. 2. — The profile of tlie central source at 2.2 /xm and 1.6 //m compared to a TINYTIM 
generated point source. 
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Fig. 3. — The plot shows the fitted spectrum (F^ oc z/^^'^) of the central source along with 
the near and mid infrared data. These assume foreground extinction by a screen so both the 
near and mid infrared points have been extinction corrected. The visible light polarimetry 
measurement (triangle) of W95 is shown as a lower limit because of the large uncertainty in 
the value. The asterisk represents the extinction corrected RL78 measurement. 
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Fig. 4. — Spectrum of an ADAF fit to the data on NGC4258. There are three components 
of emission, a Lj, oc z/^'^ synchrotron component, a L,^ ex v^^'^ power law component from 
Compton upscatter of the synchrotron photons and a bremsstrahlung component with an 
exponential cutoff at frequencies hz/ > kT^. The components from v > 10^^ Hz agree with 
M97. Below that we have adopted a z/^^^ synchrotron component instead of the z/^/^ predicted 
in M97. The W95 visible light polarimetry is shown as a lower limit because of the large 
uncertainty in the measurement. The plot also shows the 1.1 fim. upper limit derived from 
our observations, the 22 GHz upper limit of Herrnstein et al. (1998) and the Makishima et 
al. (1994) X-ray measurement. 
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Fig. 5. — Color correction between NICMOS magnitudes and ground-based magnitudes 
based on observations of standards of different colors. The lines are the best fit first-order 
polynomials described in the text. The color correction at [2.2] is small enough that we take 
it to be 0. 
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Table 1. Energy Distribution for the Compact Source 



A 


u 


Observed F,^ 


Extinction Corrected Fj, 


Reference 






Hz 


mJy 


mJy 






2-10 keV 


1.45x10^^ 




4.7x10-^ 


Makishima et al. 


(1994) 


0.55/im 


5.45x10^^ 




1.7-170 


W95 




1.1/im 


2.73x10^^ 


<0.5 


10 






l.Gfira 


1.88x10^^ 


0.9±0.3 


16 






2.2/im 


1.36x10^^ 


4.0±0.7 


25 






3.45/im 


8.70x10^3 


20±3 


46 






lO.Syum^ 


2.86x1013 


100±12 


195 


RL78 




12.5/im 


2.40x10^3 


165±20 


270 






17.9;um 


1.68x10^=^ 


300±30 


435 






1.4 cm 


22x10^ 




<0.22 


Herrnstein et al. 


(1998) 


^8-13 fim 


. filter. 











Table 2. Photometry of Unextincted NICMOS Stars*' 



Object 



[1.1] 



J 



[1.6] 



H 



[2.2] 



K 



G191-B2B (WD) 12.49±0.02 12.68±0.02 12.71±0.02 12.72±0.03 12.83±0.02 
P330E (G star) 12.01±0.02 11.96±0.02 11.57±0.02 11.59±0.02 11.49±0.02 11.51±0.02 
BRI0021 (red) 12.45±0.02 11.84±0.01 11.28±0.02 11.09±0.01 10.48±0.02 10.55±0.01 



'Obtained from the STScI wesite with the errors not representing calibration errors. 



